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accuracy, the stockout cost element displays a 
hypersensitive reaction. 
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ABSTRACT 


An analysis is made of the Navy's demand forecasting 
process and its impact on inventory system effectiveness. 
The current Navy Uniform Inventory Control Point (UICP) 
forecasting model is compared with an alternative computer- 
oriented technique using UICP data. The comparison high- 
lights the presence of highly erratic patterns in the UICP 
demand data base. Next, a simulation model is exercised 
to suggest how the UICP demand reporting method might 
contribute to the variance of recorded demand. The thesis 
concludes with another simulation indicating the relation 
Seeaemand forecasting accuracy on each component of total 
inventory cost. This simulation suggests that, while holding 
and ordering costs remain relatively insensitive to fluctuea- 
tions in forecast accuracy, the stockout cost element 


displays a hypersensitive reaction. 
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i. INTRODUCTION 


imventeory GControl 1S agpivotal activity of any logistics 
Organization. Multi-item inventory systems encompass trade- 
offs in balancing customer service needs with operating 
costs. This management task is particularly challenging 
in the military setting where item availability often 
affects mission readiness. Clearly, one objective of any 
inventory doctrine is to succeed in making those decisions 
which minimize operating costs while providing an acceptable 
level of service for a forecasted rate of demand. The 
demand forecasting process and its influence on inventory 
system effectiveness are the subjects of this study. 

First, some popular forecasting methods are profiled in- 
terms of six key evaluation measures. 

In recent decades, a wide variety of forecasting methods 
has emerged. Generally, they can be assigned to one of two 
taxonomies: qualitative or quantitative. Qualitative 
techniques are regarded as the more subjective of forecasting 
approaches. Usually conducted in a setting where historical 
data is unavailable, this class often employs expert opinion 
in constructing a forecast. The Navy Supply System uses 
such a method, called Best Replacement Factors, when esti- 
Mating initial stock levels for a new item of inventory. 


In contrast, quantitative methods make extensive use of 





historical data. There the data serves as input for various 
types of mathematical models which compute the required 
forecast. Not surprisingly, advances in computer technology 
tend to popularize the quantitative-oriented forecasting 
methods. Rather than spending considerable time synopsizing 
the more common forecasting techniques, the interested 
reader is invited to consult the existing literature for 
Supporting detail (see Makridakis and Wheelwright [Ref. 1]). 
However, it is interesting to classify a few of the widely 
igea Eechnigwues in terms of their cost, accuracy, type, 
applicability, data pattern and time horizon characteristics. 
Adapted from Wheelwright and Makridakis [Ref. 2], the pro- 
file provided by Tables 1 and 2 assists in selecting the 
most appropriate method for a given forecasting requirement. 
For example, the widespread appeal, among industry and the 
military, for the enone at smoothing technique is 
apparent. To utilize exponential smoothing a manager need 
have only three data elements: the most recent observation, 
the most recent forecast and a weighting parameter. This 
data storage consideration has been of primary importance 

to multi-item inventory systems where demand forecasts are 
routinely prepared for several thousand items. Combining 
such features as low data processing and storage costs to- 
gether with high applicability, exponential smoothing appears 
as a rational choice for the Navy's forecasting method. In 


contrast, the Box-Jenkins technique is not suitable for 
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UICP use. Although it offers great accuracy, the Box- 
Jenkins method's excessive data procesSing and data storage 
requirements render its costs prohibitive. 

However profiled in terms of these six factors, the 
best measure of effectiveness for a forecasting method is 
the economic benefit it provides to its dependent inventory 
control system. The identification of a method which incurs 
low data processing costs while contributing to minimal holding 
and stockout costs remains the goal of the inventory fore- 
caster. This notion will be pursued further. 

Before examining the cost effects of forecast accuracy 
several issues will be analyzed. The following chapter re- 
views the forecasting model currently installed at one of 
two Navy Inventory Control Points, the Navy Ships Parts 
Control Center (SPCC), and also introduces an alternate 
computer-oriented model called, "Focus Forecasting.” Chapter 
III reveals that items having erratic demand patterns con- 
stitute an appreciable portion of the ICP's inventory 
Pepulation. Chapters IV and V identify both the cost 
effects and some sources of extreme variability in the UICP 
demands. Lastly, Chapter VI uses a computer simulation to 


examine the cost effects of forecast accuracy. 
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It. TWO FORECASTING MODELS 


A. EXPONENTIAL SMOOTHING 

As previously indicated, exponential smoothing satis- 
fies the requirements for a forecasting procedure which 
demands low computer storage and run time. Accordingly, 
the UICP model utilizes a refined version of exponential 
smoothing in establishing, among other things, the parameters 
for the leadtime demand distributions. Employing a moar 
cation of the continuous-review inventory formulation found 
in Hadley and Whitin [Ref. 3], the UICP model applies the 
forecasted leadtime demand in determining both economic 
order quantities and reorder points. 

Exponential smoothing methods, originally advocated by 
Brown [Ref. 4], are a geometrically weighted sum of all past 
demands with the greatest weight applied to the most recent 


observation. Mathematically a forecast is calculated as: 


NEW FORECAST = LAST FORECAST + WEIGHT (LAST OBSERVATION 


-~ LAST FORECAST) 


Notationally this may be stated as 


which reduces to 


is 





ee te ee ee = 2 le, 
Here a@ is known as the smoothing weight and is normally 
assigned a value between zero and one. The accuracy of the 
exponentially smoothed forecast depends strongly on the 
chosen a value. In normal practice, the a smoothing con- 
stant is either selected arbitrarily or suggested through 
exhaustive sensitivity analyses. With D, again denoting 
the actual demand observation recorded in period t, the 
above equation can be represented in recursive form as: 

3 


+ a(l-a)~D eee 


+ a(1-a)“D,_ - 


F = aD, + a(l-a)D,_ 


t+1l 2 1 v2 


From this recursive form it 1s easily seen how exponential 
smoothing dilutes the effect of the older observations. 

The data-processing simplicity of this technique is evident 
Since one stored value, Fea replaces the entire block of 
t demands. One point of concern can be the initial or seed 
value for F_. Since a F, value is needed before the next 


E ie 


F forecast is prepared, a seed value must be found to 


t+1 
initialize the process when t = 0. Several solutions to the 
seed selection problem are available. Some of the simpler 
solutions recommend using the first observation as a seed or 
dividing the time series data into two parts with the first 


part reserved for initial estimation purposes (average, 


least squares estimation, etc.). In practice the seed problem 
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is only of theoretical concern. (See Makridakis and 
Wheelwright [Ref. 1].) Generally, the forecasting process 
will be in operation long enough to suppress any dependency 
@n the seed. 

The simple form of exponential smoothing described above 
is essentially an approximation of a moving average fore- 
casting process. An inherent drawback of this model is that 
it is relatively insensitive to recent trend changes. This 
weakness, left uncorrected, would frequently result in 
biased forecasts. Recognizing this limitation, the UICP 
model refines the process by incorporating two types of 
demand filters. (See Basic Inventory Manager's Manual [Ref. 
5].) First a trend test is conducted to detect sustained 
eae) Changes frem the wnderlying pattern. The trend 


statistic, TR, is a ratio test consisting of: 


2(SUM OF LAST TWO OBSERVATIONS) 


aR (SUM OF LAST FOUR OBSERVATIONS) 


A trend is considered present when either: 
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When either of these trend conditions exist, the smoothing 
constant is modified from its usual .10 value to a new, 
"heavier" weight of .30. This causes the next forecast to 
be more directly influenced by the most recent observation. 
Next, a second filter is used to check for outlier observa- 
tions. This filter computes a control tolerance band around 
Fy using multiples of mean absolute deviation (MAD), the 
average of the absolute difference between actual and fore- 


casted demand. For further discussion on the exact UICP 


use of MAD see the Basic Inventory Manager's Manual [Ref. 


5]. If the most recent demand observation lies within a 
band of width 7.5 MAD about Fie the process 1s considered in 
control and no modifications are necessary. If a single 


eut-of-control condition has been indicated by the filter, 
the outlier demand is ignored and the forecast is left 


unchanged, i.e., F = F Further, if two successive 


t+1 ol 
demand observations lie on the same "side" (high or low) 
of the tolerance band, the next forecast is calculated as 
the average of the two cohort outlier demands, l.e., 

a - Di). This condition is known as a "step 


increase/decrease." 


134 7D 


t+1 7 


This completes the discussion of the SPCC exponential 


smoothing model. Next, "Focus Forecasting" is introduced. 


B. FOCUS FORECASTING 
Focus forecasting is a new forecasting approach first 


advocated by Smith [Ref. 6]. His concept requires the 
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dynamic simulation capabilities of modern computers in 
preparing each forecast. The methodology is somewhat 
inviting due to its overriding simplicity. Focus fore- 
casting emphasizes a straightforward, flexible design in 
acquiring user understanding and confidence. The require- 
ment for transparent forecasts which are derived from simple 
Strategies and which capitalize on recent advances in com- 
puter technology motivated Smith in creating the concept. 
The mechanics of the process consist of four operations: 
backforecasting, selection, application and repetition. 
First, employing a dynamic evaluation routine, the computer 
identifies from a corpus of simple strategies the one which. 
would have best forecasted the preceding period's demand. 
Next, this selected strategy is used in preparing the up- 
coming period's demand forecast. Lastly, the process is 
repeated until "optimal" strategies have been identified and 
applied for each inventory item. Recent strides in computer 
technology make possible Smith's procedure which, when imple- 
mented for large inventory systems, requires great internal 
processing speed. For clarity, a conceptual example of the 


focus forecasting algorithm is presented. 


- Item Demand History: 
Time Period (Quarters): a Z 3 4 5 6 7 8 
Units Demanded (D,) : Buomeoce. 527 259» 270 363 357 250 


Required: Fy where t = 9. 
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Routine: 


Step 1l-~Backforecast for t-l, 1.e., compute F 


t-1 
Sera) F, = De 
= 357 
STRAT(2) = Fg, = (D_xD,)/D, 


= (357 x 259) /527 


= 175 


Step 2--Select item--"“optimal" strategy 


DIFF(1) = ABS(D, = oka, (i) ) 
= ABS({ 250 — 357) 
= 107 

DIFF(2) = ABS (D, Sl aeZ)) 


= ABS(250-175) 
= 75 


Ssanee 7/5 <e8107 item “optimal” strategy is STRAT(2). 


Step 3--Apply selected strategy 


Ee = (D 


é Cee a es 
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ry 
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9 ~ Dg XD5/D, 


(250 270) 7250 


= 261 


Step 4--Repeat for next item 


{go to step l) 


Smith has implemented a seven-strategy version for a major 
commercial wholesaler of hardware. There the system pre- 
pares demand forecasts for over 100,000 items each month. 
Smith maintains that his procedure of adapting a series of 
forecasting approaches to item demand will significantly 
Sueperform a single-formula process such as the UICP exponen- 
tial smoothing model. Before this claim is examined, we 
briefly outline how focus forecasting might be successfully 
applied to military use--specifically leadtime demand 
parameterization. 

As an initial, albeit crude, attempt at implementation, 
six simple strategies derived from three general forecasting 
categories are to be used. The first two strategies come 
mom the so-called, "Naive" class. These are the previously 
described "last period this year" method and its cohort, 
"this period last year." (Notationally--STRAT(1) = Dy! 
SMeRAT (2) = D_4:) The first method acknowledges trend; the 
second seasonality. The next three strategies are selected 
from the class of moving averages, all of which model hori- 


zontal demand patterns. The moving averages will be computed 
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as 2, 4 and 8 period averages. (Notationally--STRAT(3) = 


pe(D,_; + D,_5): STRAT(4) = .25(D,_, + -.- + D,_4); 


STRAT(5) = .125(D + ... + D, .).) Finally, as proposed 


t-39 


by Bates and Granger [Ref. 7], a composite forecast is 


t-1 


adapted in order to improve forecast accuracy by capturing 
information from each forecast strategy. (Notationally-- 
STRAT(6) = .2(STRAT(1) + ... + STRAT(5))). 

The next chapter compares the UICP exponential smoothing 


and focus forecasting models using empirical data. 
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III. MODEL COMPARISON 


A. EMPIRICAL DATA 

To evaluate the two forecasting models under considera- 
tion, a nine year demand history was retrieved. The Opera- 
tions Analysis Department at the Navy Fleet Material Support 
Office (FMSO) followed a stratified sampling procedure in 
preparing a random sample of 522 repairable and 4530 con- 
sumable items. Each of the 5052 items was represented by 
one master data record followed by several subrecords. The 
master record contains descriptive information such as 
national stock number, replacement price, etc. Each sub- 
record contains demand quantity and demand Julian date 
information. The complete record format can be found in 
[Ref. 8]. For convenience, sequential dates ranging from 
0001 to 3285 replaced the Julian dates for the nine year 
period. A FORTRAN computer program aggregated the subrecord 
demand data into 36 quarterly "buckets." (A complete listing 
of the principal FORTRAN IV source codes used in this thesis 
1s available in Appendix A). Negative demand quantities 
were assumed to be the result of customer cancellations. 
When negative demands were encountered, the quarterly demand 
balance was appropriately reduced. When adjusting for such 
cancellations, however, the quarterly balances were forced 


to maintain strictly non-negative values. 
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me EVALUATION CRITERIA 

One aim of a forecasting process is the minimization of 
the total forecast error incurred over time. When stochas- 
tic demand rates are involved, a natural assumption is that 
superior inventory control requires very accurate forecasting. 
This assumption seems reasonable since, in the extreme case 
where demand rates are deterministic, existing inventory 
models achieve zero stockout cost and optimal ordering and 
holding costs. (See [Ref. 3].) Unfortunately, forecast 
accuracy lacks an absolute standard of measurement. In- 
stead, the decision-maker is free to choose from a wide 
variety of evaluation schemes. Conceivably, each scheme 
could rank the forecasting models differently. This section 
briefly describes the evaluation criterion selected for 
comparing the focus forecasting and SPCC exponential smoothing 
models. Chapter VI presents an alternate effectiveness 
measure--the cost impact of forecast error. 

Traditionally, squared forecast errors (SFE) have been 
a useful determinant of accuracy with mean squared error 
(MSE) serving aS a popular choice. MSE is defined mathe- 


matically as: 


MSE is functionally related to variance and also enjoys 


wide acceptance as a measure of "closeness." Squaring a 
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forecast error provides two advantages. ee the alge- 
braic sign of the error is disregarded. This avoids distor- 
tion caused by offsetting positive and negative errors. 
Second, the squaring operation penalizes large forecast 
errors. 

SFE measures are not without drawbacks. For example, 
what is the preference ordering for two models posting 
MSE's of 81.49 and 75.19? The second forecasting model 
appears more accurate but with what significance? Another 
drawback of MSE is that it does not facilitate comparison 
across different time intervals. These limitations not- 
withstanding, three SFE measures will be used: mean, standard 
deviation and the aor quantile. 

The next effectiveness criterion is designed to check 
for biased forecasts. Bias distorts forecast accuracy 
through systematic overforecasting or underforecasting. 
Several statistics drawn from the distribution of forecast 
errors help identify an unbiased forecasting model. Any 
reasonably unbiased model will exhibit two attributes. 
First for large samples, the forecast errors should have a 
mean of zero. Second, as an indication of symmetry the 
median forecast error should approximate its mean. The 
standard deviation of forecast errors will also be included 
aS a measure of dispersion. 

A third important consideration in determining forecast 


accuracy is how well the model performs for the higher 
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priced inventory. The manager may benefit little from a 
model providing low bias and dispersion for only the lower 
priced items. Therefore, forecast errors will be weighted 
by item replacement price (PWFE) and tested by the three 
bias measures described above. 

Finally, the correlation between forecasted and actual 
demands 1s considered. Often two time series are closely 
related but not in a statistically dependent sense. The 
pearelation coefficient can measure the strength of the 
linear relationship between two random variables. Also 
known as Pearson's product moment coefficient and denoted 


by r, the correlation coefficient is defined as: 


where D is mean actual demand; F is mean forecasted demand; 
mee {-1,1]. 

A forecasting mechanism which parallels changes in the 
demand pattern is highly desirable. A high r value would 
typically indicate such behavior although it must be re- 
garded with caution. For example, a high r value may 
indicate a spurious relationship caused by chance or by 
the elimination of a third explanatory variable. Nonethe- 
less, r remains a recognized measure of the tracking capa- 


bility of a forecasting model. 
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In summary, four classes of evaluation criteria will 

be reported: 

1) SFE--(MSE,SD,Po9)) as measures of closeness, 

2) FE--(mean,median,SD) as measures of bias, 

3) PWFE--(mean,median,SD) as measures of price weighted 

bias, 
4) CORRE--(r) as a measure of association and direction. 
Before reporting the results of the model comparison, 

some amplifying remarks are offered. The SPCC exponential 
smoothing model was translated into the ES FORTRAN source 
code found in Appendix A. This code contains all of the 
provisions described in Chapter II. The seed values were 
computed as the average of the first four quarters of demand. 
To dilute the effect of this selection, the program was 
stabilized over demand years two and three. The demands 
representing years four through seven were examined in terms 
of the selected effectiveness measures. The FOCUS FORTRAN 
source code represents the previously described focus fore- 
casting model. For the six simple forecasting strategies 
selected, neither initialization nor stabilization was re- 
quired. To preserve an equal footing, however, demand years 
four through seven were again used for comparison to the 
SPCC model. Note that each strategy must require no more 
than eight quarters of demand history. This constraint is 
imposed since UICP data files presently access only two 


years of demand history. Also note that the FOCUS program 


27 





lacks a high demand filter. Some basic filters are offered 
by Gmith bute none were coded for the initial trial run. 

The tables which follow show the results of testing the 
SPCC exponential smoothing and focus forecasting models 
using four years of actual UICP demand history. To facili- 
tate model comparison between repairable and consumable 
inventory items, the demand history was separated into two 
populations. These tables record each forecasting model's 
performance in terms of the previously described forecasting 


effectiveness measures. 


C. FINDINGS 

The results in Tables 3 and 4 do not suggest a clear 
modelling preference. However, they do reveal some unexpected 
results which led to a redirection of this study. First 


the MSE and SFE P values indicate that the presence of 


90 


large outliers severely distorts MSE. The use of Poo as 

a more stable measure does not indicate a significant dif- 
ference in model closeness. Second, both models succeed 

in providing unbiased forecasts with the focus forecasting 
model generally producing a smaller price weighted forecast 
error. Next, there does not appear to be any significant 
difference in the model's tracking ability as measured by 


r. However, both models appear more successful at tracking 


consumables than repairables. 
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The next originally planned task was to fine tune focus 
forecasting using high demand filters and inject alternate 
forecasting strategies for subsequent trial runs. However, 
this task was interrupted when demand patterns, such as the 


three actual demand histories shown below, were observed. 


DEMAND HISTORY 


QUARTER 
ITEM IE Z 3 ~ 2 6 7 8 
1 76 5843 18798 iS 58 746 no 0 
Z 6 2 6 78 0 8 LO 3 
3 Jus 44 64 7 100 604 130 2239 


It was startling to observe a significant amount of varia- 
bility in both the repairable and consumable demand histories. 
To gain a rough estimate of the magnitude of this dispersion 
the average coefficient of variation (STANDARD DEVIATION [D]/ 
EXPECTED VALUE [D]) was computed. The two data sets 


recorded average coefficients of: 


Repairable - COEFF VAR = 3.80 
Consumable - COEFF VAR = £3.58 


Since coefficients of variation greater than 1.0 are usually 
regarded as an indication of excessive dispersion [Ref. l], 
it appears that the demand histories profiled above are the 
rule rather than the exception. Recognition of the data's 


extreme variability is further suggested by SPCC's 1980 


Ze 





change to their model's high demand filter. The control 
tolerance band was revised from its previous Fy + 3.75 MAD 
to a widened range of FL + 7.50 MAD. These erratic demand 
patterns challenge the foundation of those inventory models 
(such as SPCC‘'s) which assume a stationary demand rate from 
quarter to quarter. Thus, two new research directions were 
suggested: (1) identify the causes of erratic UICP demand 


and (2) assess the cost impact of forecast inaccuracy on the 


type of inventory models implemented by the ICP's. 
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IV. THE ECONOMIC BENEFITS FROM VARIANCE REDUCTION 


Before studying the causes of erratic UICP demands, it 
iS proper to first identify the marginal benefits which can 
be derived from reducing demand variance. 

The stochastic backorders model developed in [Ref. 3] 
contains three cost elements: order, holding and stockout 
costs. Mathematically these are represented in the total 


annual varieble cost equation as: 


co 


K(Q,r) = AA/Q + IC(Q/2 + r-u) + TA/QE f xh(x)dx - rH(r)], 
hit 


or 


hye = WerbER teHOLDING + STOCCKOUT COSTS 


For our purposes, K(Q,r) is assumed to be both differen- 
tiable and jointly convex in Q and r. Under these assumptions, 
the values Q* and r* which minimize K(Q,r) are determined 
by the methods described below. It should be noted, however, 
that joint convexity depends on the particular distribution 
of h(x). Brooks and Lu [Ref. 9] and Veinott [Ref. 10] 
show that K(Q,r) can be nonconvex when h(x) is non-decreasing 
and r < uw. (Nonconvexity may lead to a failure in the 


Speimization technigue which follows.) 


oes 





Under joint convexity the optimal values Q*, r* (where 


ie OF °< o>, 0 < r* < ©) must satisfy the equations: 


OK _ eve = Xi = aa 
30 = OR = axtlA.+ Meee is i 
+ Q(r) = [2A(A¢na(r)/tc]?/? 
aK _ Bewley ¢ fee 
ss = 0 F a az! J Ge) inode) ] + EC 
= Te [-H(r) ] + IC (using Leibniz's Rule) 
eC (0) ee Se OC TX 


(note tm must be sufficiently large such that H(r*(Q)) < 1). 
Computing a numerical solution for QO*(r) and r*(Q) poses a 
problem of composite dependency. That is to solve for Q* 
we need to know n(r) and thus r. Secondly, to solve for r* 
me-aquires a Knowledge of QO. As a practical and accurate 


alternative, a five-step numerical iteration routine is 


invoked. 


meecheck £O6r a unique solution, 1.e., if 
ys ; ze 
[2A (Attu) JIC] tween URI gue ©* ,r 


2)) Assume n(r) = 0 and starting with i = 1, 
compute Q. = POV Waele 


Be 





3) Compute rs - rem H(r,(Q.)) = Q,IC/Tr by “consulting 
(tee di Stitoueron fables of H(x) 
4) Compute OF) by using rs EQumchiind n(r,), Liew 
2 = i a2 
compute O54) = [2A (Atmn(r,))/IC] 


converge on Q. 
5) Repeat steps 3) and 4) until Oss G41 g Qi,r,. 


. Z 
To determine the effect of demand variance (0 ) on K 


requires a knowledge of dK/30 where 
ere, 6) = AA/O(r,0) + IC(O(r,c)/2 + r(o)-y) + mAn(r,9)/Q(r,0) 


Unfortunately, developing a determinant form of dK/do is 
more complicated than the earlier task of solving for Q* and 
r*, The complication stems from the composite dependency 
among go, Q and r. That is, 3K/3o0 depends on 3n(r)/30 which 
depends on r which in turn depends on Q. For example, if 
the leadtime demand distribution is assumed normal with 


mean wu and standard deviation o then 


on (e) 
a0 





sal J (x-r) 6 (x,0) ax] 


“J = = o(x,0)dx + el hoee) ae do 


where $(X,0) = (1//2mo) exp { (x-u) “/207] and r is the solution 
@emH{r,0) = ICQO(r,0)/TA. 
To remedy this complication and to graphically illustrate 


the relation between K and o, a sensitivity analysis was 


oe 





conducted using the FORTRAN program entitled INVENTOR. For 
a set of input parameters which remained fixed throughout 
the analysis 1000 realizationsof o were used in computing 
the three cost elements of K. The particular set of input 
parameters (A = 500, A = 180, IC = 15, mw = 500, uw = 90) 
was selected to ensure joint convexity in Q andr. The 
1000 values of o are equally spaced over the interval [10,350]. 
The upper and lower limits of this interval were selected to 
reflect the coefficients of variation found in the UICP 
sample. The upper limit of this range was specified to 
coincide with the average coefficient of variation (approxi- 
mately 3.80) reported in Chapter III (i.e., 350 = 3.80 x90). 
Graphical representations of the sensitivity analysis are 
found in Figure 1-4. To Ree omeen ect tson of the changes 
in each of the four cost categories, each axis shown in 
Figures 1-4 was normalized in the following manner: 
(1) Ordering, holding, stockout and total variable costs 
were each computed 1000 times for the 1000 different 
values of o. 
(2) The first value of each cost (where o = 10) was used 
as a basis value. 
(3) Each of the four sets of the 1000 cost values was 
then divided by its respective basis to produce the 
vertical coordinates used in the graphs. The 1000 oa 
values were normalized by dividing each by o= 10 
thereby producing the horizontal coordinates found 


in the graphs. 
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Figure 1. Ordering Cost Sensitivity to Sigma 





NORMALIZED CHANGES IN HOLDING COSTS 
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Figure 2. Holding Cost Sensitivity to Sigma 
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NORMALIZED CHANGES IN STOCKOUT COSTS 


Figure 3. Stockout Cost Sensitivity to Sigma 
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NORMALIZED CHANGES IN TOTAL VARIABLE COSTS 
i 


Figure 4. Total Cost Sensitivity to Sigma 
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These graphs highlight, in a general sense, the theoreti- 
cal cost effects of variability in leadtime demand. Figure 
1 shows that as o increases over the relevant range the 
ordering costs developed in the Hadley-Whitin <Q,r> model 
decay in a nonlinear manner. In contrast, Figures 2 and 3 
show the corresponding effects on the holding cost and 
stockout cost components of the <Q,r> model. There both 
cost elements increase in a roughly linear fashion over the 
relevant range with (d3BC/d0) > kx (dHC/d0) (where k = 4.5). 
The reaction of total variable costs to o has been included 
to show the overall effect. It should be remembered that 
the effects of each of the other three cost elements will 
be dampened when translated into total variable cost. 

This sensitivity analysis indicates that a study of the 
sources of leadtime demand variance would help reduce holding 
costs but would also generate particularly beneficial effects 
on stockout costs. The identification of the causes of 


demand variance is the subject of the following chapter. 
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V. Shy eeriesezeD CAUSES OF DEMAND VARIABILITY 


Demand randemness results from two sources. There is 
the truly patternless variation that is unpredictable and 
uncontrollable, but there is also the variation which re- 
wares trom the way data are collected or reported during a 
given operational schedule. This second source of variation 
is the subject of analysis of this chapter. 

Erratic items constitute an appreciable portion of the 
Military's inventory population. Silver [Ref. 11] addressed 
this issue and it is also demonstrated by the coefficients 
of variation found in the UICP sample data discussed above. 
Silver maintains that imperfections in the wholesale manage- 
ment information reporting system are a primary cause of 
Peak ic demand. (FMSO calls their inventory reporting network 
the Transaction Item Reporting (TIR) System.) 

The TIR system is a tri-level organization which sub- 
scribes to a policy of geographical responsibility for supply 
Support. The first level of the TIR's multi-echelon network 
is the end-user or consumer level (e.g., a shipboard supply 
department). Level two is composed of several retail out- 
lets such as the at-sea Mobile Logistics Support Force 
(MLSF) ships and the shore-based stock points (Naval Supply 
Centers, industrial rework facilities, Naval Air Stations, 
etc.). The final tier contains the Navy's two wholesale 


activities--the Inventory Control Points. 
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Under the most common scenario, a shipboard supply 
department receives a routine demand request from one of 
its customer departments. If the demand reduces the ship's 
available stock past the reorder point or if the demand 
is for an unsupported item, a replenishment requisition is 
submitted to the nearest retail outlet by one of two means. 
For ships operating at sea with an assigned MLSF unit, 
the replenishment requisition is submitted directly to the 
MESF ship. For ships operating without MLSF support or for 
those which are inport when the demand occurs, the replenish- 
ment requisition is deposited with the nearest shore-based 
retail outlet. Clearly for situations where a customer-ship 
operates at sea for extended periods without MLSF support, 
the replenishment requisition encounters substantial delays 
before reaching the TIR's retail level. Under either case 
when the requisition ultimately reaches a retail activity, 
the TIR system activates and provides transaction status to 
all levels. When the available stock at the retail outlet 
reaches its reorder point, a collective replenishment requi- 
sition is submitted to the ICP. This requisition reflects 
the current and anticipated future demands from the retail 
outlet's geographically assigned customers. In this cas- 
cading fashion every demand which originates at the consumer 
level ultimately reaches the wholesale level where the ICP 
is tasked with the system-wide reorder and budgetary 


responsibilities. 
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A Simulation was conducted to illustrate how the TIR 
Organizational structure may cause truly non-erratic con- 
sumer level demands to be observed as the highly erratic 
demand actually seen at the ICP level. (See Lewis and Uribe 
[Ref. 12] and FORTRAN source code, DEMSIM.) The framework 
for this simulation was adopted from Shields [Ref. 13]. 

The simulation modified the present TIR demand reporting 
scheme by inserting an artificial reporting channel which 
reports all consumer demands directly to the ICP as they 
eur. To date this reporting medium has been reserved for 
only high priority requisitions. This alternate reporting 
channel provides two valuable capabilities. First by directly 
reporting user-demands, the ICP‘s records will accurately 
reflect the demand patterns occurring on the consumer level. 
Second, the ICP recording operation is expedited by elimin- 
ating the srt ias delays incurred when ships operate with- 
Out MLSF services. Together these accuracy and timeliness 
considerations can have a big impact on the effectiveness 

SE an ICP. 

The structure for the DEMSIM program is as follows: 

1) assume a 30 ship population with each ship operating 
independently over an18 year active service life. 

2) initial ship inport and at-sea assignments were 
randomly selected. (The terms "inport" and "at-sea" 
are meant only to describe when a retail supply 


activity is readily available). 
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3) assume inport periods follow a Uniform [5,30] dis- 
tribution and at-sea periods follow a Uniform [5,45] 
GCrUsSi=erDubaon. 

4) generate demands in a random manner following a 
compound Poisson process where demand interarrival 
times are Exponential and demand requisition sizes 
are Poisson. 

5) record demands according to: 

a) the indirect reporting system presently used 
by TGP. 

b) the artificial system where the ICP receives 
observations directly from the consumer. 

6) aggregate each demand record into quarterly groupings 
and compute the variance of each time series. 

7) repeat the process 50 times. 

The result of the simulation supports a similar finding 
to that of Shields [Ref. 13], that demand variance is sub- 
stantially reduced using the alternate reporting method. 
After 50 replications the direct reporting scheme demon- 
strated an average variance of 3050 as opposed to 5560 for 
the indirect ICP reporting system (a 45% reduction). 

This simulation is meant to highlight only one cause of 
erratic demand. Presumably other changes in the demand 
recording system could yield additional improvements. 
Examples of such changes are: 

* Limited Demand History~-continued strides in computer 


technology will make storage of longer demand histories 
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an economic reality. Greater depth in the demand 
history would facilitate aggregation by other than a 
quarterly basis. A change in the base time period 
could provide a smoother stream of observations. 
Customer Identity--The interaction of large customers 
(e.g., naval bases) who order more product less often 
with small customers (e.g., destroyers) who order 

less product more often might cause high variability. 
It may prove advisable to monitor a few large 

customers on an individual basis. 

Customer Buying Habits--A practice of batch ordering 
could easily be a contributing factor. Actual consump- 
tion of material may follow a stable pattern, however, 
due to such factors as funding and increased operational 
commitments material is often ordered in large lots. 
Product Identity~--The disaggregation of products 

into their life cycle states (growth, steady state, 
gradual phase-out) might produce different and dis- 
tinct demand patterns. 


Whether such changes could or even should be implemented 


largely depends on the results gleaned from comparing the 


inventory cost savings against the cost of altering the 


demand recording system. This study does not undertake 


such a task but focuses instead on a more pertinent issue-- 


the dollar cost of forecast accuracy. That is, for a given 


demand pattern how desirable is it from an inventory cost 


perspective to pursue a more accurate forecasting model? 
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Vi. COs eacisor FORECAST ERROR 


The final phase of the study uses a computer simulation 
to examine how forecast error affects ordering, holding and 
stockout costs. This examination does not suggest how to 
improve forecasts but it does indicate if it is worth doing 
anything at all. 

The <Q,r> model presented in Chapter IV was a repre- 
sentation of a theoretical inventory cycle for which sta- 
tionary annual demand rates were applicable. Here the 
interest is in meaSuring the actual inventory costs which 
an ICP incurs when annual demand, due to forecasting 
inaccuracy, is no longer stationary and when both the 
reorder quantities (Q;) and reorder points (r.) vary for 
each oe cycle. To conduct such an analysis requires a 
type of "bookkeeping" computer routine capable of both 
maintaining the three inventory cost "accounts" and com- 
Pucing Q; and rc. for each cycle. 


es cycle length 


For the purposes of this chapter, the i 
(T.) is defined as the interval of time between receipt of 
successive orders. Registering cycle order costs reduces 
to the simple task of posting the cost per order (A) once 
each cycle. The second account, holding costs, is posted 
in one of two ways. First, if no backorders occur then the 


cycle holding costs are ST. [ (Opening Inventory) ; + (Closing 


Inventory). ,]. If a backorder does occur, however, the 
i-l 
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cycle holding costs are just 5(T! (Opening Inventory) | | where 

T* is that increment of T for which the inventory balance 

was poSitive. It should be noted that in exercising this 
actual inventory cycle the opening cycle inventory may 

differ from the Qt+(r-u) quantity found in the theoretical inven- 
tory cycle of Chapter IV. A more accurate representation is 
nQ,+(rj-u,) where n iS an integer. This is necessary to ensure 
that the reorder point for each successive cycle will be 
reached. That is, since Q. and rs change with each cycle 

under certain Situations an integer meneeuD ile of Q. must be 
ordered during the ae cycle to enable the i+l opening 

cycle inventory to exceed its reorder point, Ceaye 

As to the third account, stockout costs, two possibili- 
ties must be considered. The first possibility is that the 
backorder cost 1S 7 per actual backorder with no penalty 
assigned for the length of time the backorder exists. | 
Second, the actual backorder cost, T, could also be con- 
sidered as a function of the time the backorder exists 
(T, -Ti). Separate accounts are used to examine the effects 
of forecast error on each possibility. 

Before presenting a summary of the structure of the 
computer simulation, a measure of forecast error 1S required 
to demonstrate the cost effects for various magnitudes of 
forecast error. The difficulties of dealing with a single 
criterion of forecast accuracy was addressed earlier. In 


Spite of the fact that it can at times be an unreliable 
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figure of merit (see Tables 3 and 4), MSE was selected as 
the measure of accuracy. Specifically, an array of eleven 
equally spaced values ranging from 0 to 5000 was chosen 
miee., MSE(K) € 70,5000], k = 1,11). This range of values 
partially reflects the MSE quantities found in Tables 3 
and 4. 

The simulation computer program (see MSESIM FORTRAN) 
borrowed the input parameters found in the INVENTOR source 
code discussed in Chapter IV (with o = 10 here). In the 
MSESIM routine a demand value (A; ~ NORMAL(180,225)) anda 
leadtime (LT. ~ NORMAL(6(mos),1)) were generated for each 
cycle using an available random number generator (see [Ref. 
12)). Using the synthetic demand rate, the available cycle 
stock level is depleted until the cycle reorder point (r.) 
is reached. This triggers the placement of an order of size 
nQ; such that this replenishment quantity will exceed Cea: 
Stocks continue to deplete at a rate of A; over the syn- 
thetic leadtime period (LT.). Upon reaching the end of LT, 
the cycle qT. 1s completed and the reorder quantity nQ; 
arrives. At this point the accounts for ordering, holding 
and stockout costs are updated to reflect the costs actually 


ana. are com- 


Baeurred during the cycle. Next, | 


Qa! 
puted using the basic algorithm found in Chapter IV except 
that ds is deliberately falsified to reflect the forecast 
error introduced by the chosen forecasting model. (Fore- 


cast errors (FE. ) were randomly generated as NORMAL(MSE(k), 


1) deviates.) In this way, the FE. provided by the 
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forecasting model causes a non-optimal Qs44 and ed Eo 

be computed. This process continues until a time counter 
reaches a specified upper limit at which time the evolution 
is rerun for another pass. After fifty passes the next 
candidate forecast error (MSE(K+1)) 1s introduced and the 
entire simulation process repeats for another set of fifty 
passes. After all eleven values of MSE(K) are used the 
computer program plots the grand average for each cost 
account as a function of its corresponding value of MSE(K) 
producing the graphs which follow. 

The figures presented below are the results of the MSESIM 
computer routine just described. Again each ordinate axis 
is normalized to facilitate comparison. Also for complete- 
ness both forms of stockout costs (unit and time weighted) 
and their corresponding total variable cost graphs follow 
the ordering and holding cost plots. 

A review of Figures 5-10 shows that both ordering and 
holding costs exhibit relatively minor changes over the 
relevant range (9% decrease; 6% increase, respectively). 
The roughly linear increases in both stockout costs are 
guite dramatic especially for time-weighted backorders where 
greater than a five fold increase is observed. Figures 8 
and 10, the two total variable cost graphs, must be viewed 
with caution for, once again, their behavior is directly 
influenced by the particular values selected for A, IC, 7, 


O. 
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Figure 5. Ordering Cost Sensitivity to MSE 


49 





NORMALIZED AVERAGE HOLDING COST 
te 
& 


1 
@ 800 1000 $00 2000 2800 3000 3800 4000 4500 5000 
MEAN SQUARE ERROR 


Picune some Holding Cost Sensitivity to MSE 


30 





See eee: 


& 


& 


NORMALIZED AVERAGE STOCKOUT COST—UNIT WEIGHTED 


Figure 7. Stockout Cost (Unit Weighted) Sensitivity 
to MSE 
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NORMAUZED AVERAGE TOTAL VARIABLE COST—UNIT WEIGHTED 


Figure 8. 
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Total Cost (Unit Weighted) Sensitivity to MSE 
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Figure 9. Stockout Cost (Time Weighted) Sensitivity 
tO MSE 
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Figure 10. Total Cost (Time Weighted) Sensitivity 
to MSE 
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- The conclusions gleaned from these graphical displays 
are two. First, the consequences of even severe forecast 
error on both ordering and holding costs appear relatively 
insignificant. The second and more profound finding is 
that a particular forecasting model's accuracy (as measured 
by MSE) appears to play a vital role in keeping total 


stockout costs reasonable. 
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Vit. CONCLUSIONS 


A. SUMMARY AND FINDINGS 

The study began with a comparison of the SPCC forecasting 
process to an unrefined form of a newer computer-oriented 
technique. The use of several effectiveness measures showed 
neither exponential smoothing nor focus forecasting per- 
formed significantly better in forecasting actual SPCC 
demands. One outgrowth of the comparison was the prepon- 
derance of erratic-demand patterns in the ICP inventory 
population. Next a sensitivity analysis was designed to 
measure the cost effects of demand variability. This analy- 
Sis indicated that reductions in demand variability gener- 
ate generous cost savings, particularly in the form of 
reduced stockout costs. Not content with viewing the erratic 
demand situation as inflexible, two research directions were 
pursued: (1) a change to the demand pattern; (2) the cost 
effects of changing the forecasting model. First the in- 
clusion of an artificial demand reporting channel produces 
considerable reductions in demand variance. A second and, 
presumably, more complicated improvement to the demand 
forecasting process would be the development of a model 
capable of accurately forecasting erratic items. To gain 
insight into the expected economic worth of any such model, 


the study concluded with an examination of the cost effects 
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Gemtonecaset error. There it was found that stockout costs 


decrease drastically as forecast accuracy improves. 


Ee 


RECOMMENDATIONS FOR FURTHER STUDY 


Follow-on research in the areas listed below may aug- 


ment the research conducted in this study. 


1) 


3) 


Conduct further investigation into modifying the TIR 
reporting system by routing ‘all consumer demands 
dimectly to the ICP. 

Investigate other causes of demand variability and 
examine the corrective procedures necessary to 

smooth input demand data. (Several possible sources 
were hypothesized in Chapter V.) 

Continue to analyze computer-oriented forecasting 
techniques. Since the development of the exponential 
smoothing concept, technological advancements in the 
computer field have lessened computational and data 
storage costs. Newer, innovative forecasting methods 
capitalize on these advancements. Research into 
dependable forecasting models which are easily under- 
stood, easily maintained and easily adapted by a multi- 
item inventory activity should receive continued 
sponsorship. The focus forecasting concept is one 


possible candidate. 
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